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Abstract. We calculate the bulk and boundary parts of the free energy for an 
open spin-1/2 XJfZ-chain in the critical regime by bosonisation. We identify the 
cutoff independent contributions and determine their amplitudes by comparing 
with Bethe ansatz calculations at zero temperature T. For the bulk part of the free 
energy we find agreement with Lukyanov's result [Nucl. Phys. B 522, 533 (1998)]. 
In the boundary part we obtain a cutoff independent term which is linear in T 
and determines the temperature dependence of the boundary susceptibility in the 
attractive regime for T <S 1. We further show that at particular anisotropies 
where contributions from irrelevant operators with different scaling dimensions 
cross, logarithmic corrections appear. We give explicit formulas for these terms 
at those anisotropies where they are most important. We verify our results by 
comparing with extensive numerical calculations based on a numerical solution of 
the T = Bethe ansatz equations, the finite temperature Bethe ansatz equations 
in the quantum-transfer matrix formalism, and the density-matrix renormalisation 
group applied to transfer matrices. 



PACS numbers: 75.10. Jm,75.10.Pq,02.30.Ik 

1. Introduction 

The spin-1/2 XXZ-ch&in is described by the Hamihonian 

N,N-1 r. 



(1) 



Here J > is the couphng constant and N is the number of sites. The model is critical 
for anisotropy —1 < A < 1. If we take N as upper boundary in the sum in Eq. (1) we 
have periodic boundary conditions (PBC) whereas taking the sum only up to — 1 
corresponds to open boundary conditions (OBC). 

The reasons for the popularity of this model are twofold: On the one hand, 
this model does indeed capture the basic physics in some real physical systems as for 
example SrCuOa [1] . On the other hand, the model is exactly solvable by Bethe ansatz. 
Furthermore, in the continuum limit at low energies it is equivalent to a free boson 
model up to irrelevant operators. The exact solution at zero temperature T for the 
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isotropic antiferromagnetic case A = 1 and PBC has been first constructed by Bethe 
[2]. This so cahcd Bcthc ansatz (BA) has later been gcnerahscd to allow also for the 
calculation of ground-state properties in the anisotropic case [3, 4]. Thermodynamic 
properties can also be studied by using either the thermodynamic Bethe ansatz 
(TBA) [5] or the quantum transfer-matrix (QTM) approach [6, 7]. However, at finite 
temperature even the equations obtained for simple thermodynamic quantities like 
free energy, specific heat or susceptibility can often only be solved numerically. 

The BA solution for the XXZ-cha.m with OBC at T = has been constructed 
by Gaudin [8], Alcaraz et al. [9], and Sklyanin [10]. The application of the TBA 
for this case seems to be difficult or even impossible as is discussed in more detail in 
[11, 12]. A modification of the QTM approach is possible as demonstrated in [12] but 
an evaluation of the obtained formulas still seems to be a formidable task. 

An entirely different approach is based on bosonisation (see [13] and references 
therein). Here a low-energy effective theory is derived which is just a free boson model 
up to irrelevant operators. The main advantage of this approach compared to the BA 
is that any thermodynamic quantity and also any kind of correlation function can 
be calculated very easily for the free boson model. Corrections to this free boson 
approximation are obtained by doing perturbation theory in the irrelevant operators. 
A disadvantage is that the spin- wave velocity v, the Luttinger parameter K as well 
as the amplitudes of the irrelevant terms in the Hamiltonian can only be obtained 
perturbatively for |A| <C 1. However, one can use the BA to fix v and K which 
in the language of the bosonic model is equivalent to summing exactly an infinite 
series of terms which woiild renormalise the bare values of and K. Lukyanov has 
demonstrated that the BA can also be used to fix the amplitudes of the leading 
irrelevant operators in the bosonic model [14]. This allows to give analytic formulas 
at low temperatures for the free energy and derived quantities. It also allows to 
calculate the correlation amplitudes for the leading and some sub-leading terms in an 
asymptotic expansion of spin-spin correlation functions [15]. 

The XXZ-chain with open boimdarics has attracted considerable attention 
recently [16, 17, 11, 18, 12] because it is the simplest model for a spin-chain containing 
a small number of non-magnetic impurities which cut the chain into parts with 
essentially free boundaries. It has been shown that the leading irrelevant term in 
the bosonic model for anisotropy 1/2 < A < 1 leads to a contribution in the boundary 
susceptibility Xb (C(1) part of the susceptibility if the bulk part is 0{N)) which 
diverges as 7^2^-3 ^-Jjq Luttinger parameter K varies between K = 3/2 at 

A = 1/2 and = 1 at A = 1. The 1/T behaviour (up to a logarithmic correction) at 
the isotropic point is of special interest because it is a Curie-like contribution without 
any free spins being present. The isotropic ferromagnetic point has also been studied 
recently [18] and it has been found that the boundary susceptibility xb ~ —^/T^, 
i.e., the boundary susceptibility diverges more rapidly and with opposite sign than 
the bulk susceptibility x ^ 1/T^. 

In this article we want to use Lukyanov's method [14] to obtain explicit 
expressions for the free energy of the open XXZ-chain in the entire critical regime 
— 1 < A < 1. To do so we will apply perturbation theory in the leading irrelevant 
operators of the bosonic theory and compare with results obtained by Bethe ansatz. 
This will allow us to fix the amplitudes of the irrelevant terms. Our result for the 
bulk part of the free energy will confirm Lukyanov's result, however, we will obtain in 
addition the boundary part. In particular, we find that the boundary susceptibility in 
the attractive regime behaves as xb ^ const compared to the xb ^ const -|-T^^~^ 
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behaviour found before for < A < 1 [16, 17, 11]. Another focus of our work are 
logarithmic corrections. Wc will find that the amplitudes in the low-temperature 
expansion of the free energy as well as the amplitudes in the expansion of the ground- 
state energy in terms of the magnetic field h show divergences at infinitely many 
anisotropics. We will show that these divergences occur because contributions from 
terms with different scaling dimensions cross at these points. The divergences in the 
two terms which cross cancel each other leading to a logarithmic correction. We will 
give explicit formulas for these corrections at those anisotropics where they are most 
important. 

Our paper is organised as follows: In section 2 we bosonise the open XXZ- 

chain and expand the free energy in powers of h, T by using perturbation theory in 
the leading irrelevant operators. In section 3 we present the zero temperature Bethe 
ansatz solution for this model. In particular, we show how the Wiener-Hopf method 
can be used to obtain leading and next-leading terms in an expansion of the groimd- 
state energy in powers of h. We then determine the amplitudes of the irrelevant 
operators in the field theory approach in section 4 by comparing with BA. We confirm 
Lukyanov's result obtained by the same approach for the periodic case. In section 
5 we summarise our result for the bulk and boundary part of the susceptibility. We 
analyse where the amplitudes of the different terms diverge and give explicit formulas 
for the logarithmic corrections caused, at points where they appear in the leading or 
sub-leading term. We check all our analytical results by comparing with a number 
of numerical data obtained by a numerical solution of the T = BA equations, the 
BA equations in the QTM approach, and the density-matrix renormalisation group 
(DMRG) applied to transfer matrices in section 6. In the final section we give a short 
summary and conclusions. 



2. Bosonisation for the open chain 

First, we want to briefly review the bosonisation procedure for the XXZ-chain (sec for 
example [13]) to clarify our notation and to remind the reader how a priori unknown 
parameters in the bosonic model can be determined by Bethe ansatz. 

Using the Jordan- Wigner transformation we can represent the XXZ-chain with 
OBC as given in Eq. (1) in terms of fermionic operators tpj: 

H = Ho + Hini (2) 

N 

where we have added two sites and impose the boundary conditions tpo — ipN+i = 0. 
For A = we are left with a free-fermion Hamiltonian Hq. The low-temperature 
properties of this model can be studied by linearising the dispersion around the Fermi 
points zLkp and going to the continuum limit. To do so, we introduce right and left 
moving fermions by 

^ « e^'^^ijRix) + e-'^^^^L{x) (3) 

with X = ja and a being the lattice constant. In the half-filled case we have kp = /2a 
(no magnetic field applied in (1)). The continuum limit is achieved by approximating 



The open XX Z -chain 



4 



il^j+i = '4}{x)+adxtj>{x). This leads in lowest order in a to the continuum Hamiltonian 

Ho = ivo / dx tplid^ipR - tpld^ipL (4) 
Jo L ^ 

where L = Na and the spin-wave velocity is given by vq = Ja. 

This Hamiltonian can be bosonised by 

i^R = — ^ exp(iV47r</>ii) ; ipL = — ^ exp(-i\/47r(/)z,) (5) 
VZTT v27r 



where (j) = (jjR + ipL is a bosonic field obeying the standard commutation rule 
[(/>(a;), n(a;')] = \5{x — x') with 11 = VQ^dttj)- Up to total derivatives the Hamiltonian 
(4) takes the following form in terms of the bosonic fields 

Ho = vo r dx [{d^cj^Rf + {d,,6L f] = ? [tf + (a,</>)2] . (6) 

Jo ^ Jo 

The interaction part of the Hamiltonian (3) can also be expressed in terms of the 

bosonic fields 



j 

~ ^ / dx [{dM'' + {d.Hf + 'iidMidM] ■ (7) 

The first two terms are identical to the free part (6) and yields only a renormalisation 
of the velocity 

wo^t; = t;o(l + 2A/7r) . (8) 

The third term is an interaction between right and left fields. Crucially, even with 
this term included we still obtain a free boson model in terms of the fields </>, H 

H = Ho + H^nt = lJ dx [(1 - c)n2 + (1+ c)(a,</))2] (9) 

with c = 2A/(7r + 2A). We now rescale the fields by 



/AttR ' " ' " ~ ' ~ 47r V 1 - c 47r V ' n ^"^^^ 
and define the Luttinger parameter 

1 / 4A\"^/^ 4A , , 

K= ^ = 2 1 + — «2 . 11 

The bosonised version of the Hamiltonian (1) is then given by 

H=^ ["^ dx[U' + {d,cl,f] . (12) 
^ Jo 

How serious should one take the values fovv,K given in (8) and (11), respectively? 
We have derived an effective bosonic Hamiltonian for the spin chain by expanding all 
terms up to linear order in the lattice constant a. So clearly, there are higher order 
corrections to the Hamiltonian (12). In fact, there are infinitely many irrelevant 
terms which have to be added to this bosonic Hamiltonian to obtain an accurate 
description of the spin chain. Some of these terms (in fact, also infinitely many) will 
further renormalise the values of v, K. So Eqs. (8,11) are perturbative results for v, K 
expected to be valid only if A <C 1. On the other hand, the model is exactly solvable 
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by BA and we will discuss this solution in section 3. Among other things, BA allows 
to calculate the spin- wave velocity exactly. In other words, the BA allows to sum up 
exactly the infinite perturbative series for v obtained by bosonisation. In section 3 we 
will also see that the leading term in the bulk susceptibility at zero temperature is a 
function of v and K only. This allows it to determine also K exactly and leads to the 
well known results 



V = (13) 

2 arccos A 

K=^ = T. (14) 

27ri?2 TT- arccos A ^ ^ 

Note that these exact formulas agree with (8,11) in lowest order in A as expected. 
When we use the bosonised Hamiltonian with these values for v, K wc should keep 
in mind that all contributions from perturbation theory in the irrelevant operators 
leading to a renormalisation of these parameters are already accounted for. 

Finally, wc want to give the direct relations between the bosonic and the spin 
operators [13] for later use 

S- « 2^^^?^' + (-l)''const cos ^ 

Sr w exp(27riii(^) ^const cos + (-l)^const^ . (15) 

Here <j) = 4>l — (pR is the dual field and we clearly have to identify </) = </> + 2ttR and 
(p = (f)+ 1/R, i.e., (p, 4> are compact fields. 

2.1. The spin-chain in a m,agnetic field and the mode expansion for OBC 

When we add a magnetic field term Hm = hJ2f=i'^j to the Hamiltonian (1), the 
bosonic Hamiltonian (12) has to be replaced by 



H = - I dx 
2 70 



ttRv 



(16) 



according to Eq. (15). By performing a shift in the boson field 

(j)^ (i) + hx/{2i:Rv) (17) 
we can rewrite this as 

H=lj\x[Jl^ + {dAf]-L'^. (18) 

The bulk susceptibility in lowest order is therefore given by Xbuik = Kjluv. The 
magnetic field shifts the Fermi points away from kp = 1^/2 so that according to 
Eq. (3) the right- and left-moving fermions are no longer situated near these points. 
It is, however, easy to sc!c that the shift in the Bose field (17) exactly compensates for 
this shift so that the fermion fields r,^ l belonging to the shifted Bose field have 
again kp = 7r/2. The boundary conditions for the field (p in (18) are therefore exactly 
the same as in a system without magnetic field. 

We first consider the boundary conditions for the ■^-field: 

V'(O) = ip{N + 1) = 

=>V'L(o) + V'fl(o) = 

^i,L{L) + {-lfijR{L)=Q (19) 
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where L = Na and the last two relations have been derived from (3). From the 

bosonisation formula (5) and the rescaling relation (10) it follows that [(pfi{x), (f>L{x)] — 
mE? at any site x ^ Q,L, i.e., the left and right fields do not commute in general. 
This is, however, different at the boundaries where (I>r,<Pl are related according to 
(19) leading to [M^), M^)] = [ML), M^)] = and 

<^fl(O) + 0L(O) =7ri? + 27ri?n, (f)R{L) + <j)L{L) = 2TTRn' (20) 
where n is an integer, n' is a half-integer if N is even and an integer if N is odd. 
Eq. (15) allows us to fix the numbers n, n' because S^^^ = J2j ~ — (/)(0)]/27ri?. 
The total spin Sf^, should be an integer for L even and a half-integer for L odd. Fixing 
the boundary condition for (f>{0) wc therefore obtain 

0(0) = ttR ; (t){L) =TrR + 27ri?S'f^, . (21) 
This leads to the following mode expansion for OBC [19] 

0(x, t)=nR + 2nRS:^,!^ + "^^^^^^ (e-'-"^a„ + e^^^ai,) , (22) 

n=l ^ 

where a„ is a bosonic annihilation operator. 



2.2. Correlation functions for the free boson model 

From the mode expansion (22) we can directly calculate the basic correlation functions. 
As we are not interested in finite size effects but rather in the bulk and boundary 
properties of the chain in the thermodynamic limit we will consider L ^ oo. The 
results are then valid for the semi-infinite line, i.e. , a system with only one boundary. 
We start with the basic correlation function 

(0(a;i,Ti)0(a;2,r2))oBc 

= ^ (ln(zi - Z2) + ln(zi - Z2) - ln(zi - za) - ln(zi - za)) (23) 
47r 

where r represents the imaginary time. In the last line we have introduced new 
coordinates 

z = — \{x — vt) = VT — ix , z = i{x + vt) =VT-\-ix. (24) 
For PBC we would find instead 

(0(xi,ii)0(a;2,t2))pBc 

"^-=3° - ln[(a;i - X2f + - T2f\ + const (25) 
47r 

= - (ln(2;i - zi) + ln(^i - Z2)) + const . 

From these basic correlation functions one can easily obtain correlation functions as 
for example {dx<i>dx<r') by taking partial derivatives. The correlation function for the 
exponential fields can be obtained by 

- £ ^W^IS = (26) 

n— ^ ^ ^ ^ 
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and 

( exp I 1 I exp ( ±1 j ) = ( exp ( i-(0(a;i, n) ± (/)(x2, rz)) 

= exp(^--^((.^(a;i,ri)±</.(a;2,r2)f)^ (27) 

= exp (^-l^K^^XuTi)) + (02(X2,T2)) ± 2(,^(,Xi,ri)0(,X2,T2))] 

These formulas for the correlation functions will allow us to calculate the free energy 
of the model on a semi-infinite line perturbatively in subsection 2.4. 

2.3. The free fermion and the free boson model 

The free part Hq of the Hamiltonian in Eq. (2) can easily be solved by Fourier 
transform. The free energy is given by 

N 

^ = -^Ei"(i + e"'''"'") ^'^^ ^"=(ivTi)'' ^^^^ 

where we set the lattice constant a and the coupling J equal to one. The bulk and 
boundary part can be obtained by using the Euler-Maclaurin summation formula 



(29) 



F = --{N + 1) r dk\n{l + e-'5™^'=) + | (ln(l + e'^) + ln(l + e"'')) 
Jo 2 

leading to the following low-temperature expansion 

+ 1---It^- S^r' + 0{T^) + O(Te-Vr) . 
2 TT 6 360 ^ ^ ^ ' 

This means that bulk and boundary part have exactly the same temperature 

dependence if T <C 1. The same is also true for the bulk and boundary susceptibility 

[11] 

Xw.(/i = 0) = XB{h = 0) = - + Jt^ + 0{T^) . (31) 

Consider on the other hand the free boson model (12). Using the mode expansion 
(22) the Hamiltonian takes the form 

^ = ^Et + ^) • (32) 

To identify the terms which are independent of the details of the dispersion relation 
we introduce a momentum cutoff A. The free energy is then given by 



AN , AAf 

2L+^El- 



1 I ^'^^ 

l-exp(- — 



(33) 



The first part is the ground state energy Eq which is clearly cutoff dependent. We can 
rewrite the second part and obtain 

F-E "-""Pt-^") mi 
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The exponential factor in the numerator is ^ 0(exp(— A/T)) and can therefore be 
ignored if T <C 1. If we now consider the hmit N ^ oo with T fixed we obtain 

T2 ^ 1 ttT^ 
F = Eo-L—y^ — =Eo- L— . (35) 

n— 1 

This means that the free energy of the continuum free boson model does not contain 
any boundary terms. The reason for this is the continuum limit which leads to 
a[N + 1) ^ Na = L. All boimdary contributions in the field theory therefore have 
to come from higher order corrections to the free boson model. We expect that the 
same powers of T, h which are present in the boundary free energy and the boundary 
susceptibility for the free fermion model (which we will call regular terms afterwards) 
are also present for the general case — 1 < A < 1 and only the coefficients of these 
terms will change with anisotropy. By dimensional analysis it becomes clear that the 
regular terms in boundary quantities obtained by perturbation theory in the irrelevant 
operators have to involve the short distance lattice cutoff. It is therefore impossible 
to determine the coefficients of these terms within the field theory approach. We will 
come back to this point in the next section. The only boundary terms which can be 
reliably obtained from the field theory are therefore those which have no analogue 
in the bulk part. Such terms must have an amplitude which vanishes at the free 
fermion point where bulk and boundary parts show the same dependence on T, h at 
low energies. 



2.4- Perturbation theory in the irrelevant operators 



As we have already mentioned when deriving the bosonic Hamiltonian (12) from the 
original spin-model (1) there are higher order corrections. These corrections are all 
irrelevant in the whole critical regime except of the umklapp term which becomes 
marginal at A = 1. This allows it to calculate their effects systematically by simple 
perturbation theory. 

The bosonised Hamiltonian for the spin chain including the irrelevant terms with 
largest scaling dimensions is given by 



H 



dx 



+ Ai / dx cos 



a 

+ — 



dx 



ttHv 



2(f> 

'r 



(36) 



Aidx 



with a priori unknown amplitudes Ai , a, /?. Here we have separated the terms with 
integer scaling dimensions into a part where right- and left-fields are mixed and a part 
where they remain separated. This will allow us to determine the amplitudes in section 
4 by using the fact that no right and left mixing occurs at the free fermion point. The 
coefficients a{K), b{K) depend on the Luttinger parameter K and describe the relative 
weight between those terms in each bracket with different symmetries. When we shift 
the (/(-field again as in (17) the magnetic field appears in the irrelevant terms. 

We start with the umklapp term. First order perturbation theory in Ai leads to 
the correction 

2Khx^ 



f[^'> = AiT 







/ / d^X 


^cos ^ 


Jo JO 





(37) 



The open XX Z -chain 



9 



where (• • ■)o means the correlation function calculated for the free boson model. This 
is the term which has already been considered in [16] and [17] to calculate the leading 
contribution to the boundary susceptibility and specific heat for 1/2 < A < 1 at zero 
magnetic field. We will derive here the more general expression for the free energy 
contribution at finite temperature and finite magnetic field. Using (26) we obtain at 
zero temperature 

poo „„„ f 2Khx \ 

To obtain the result at small finite temperatures we can use the standard conformal 
mapping of the complex plane onto a cylinder to replace the correlation function (26) 
by its finite temperature coimtcrpart so that (38) becomes 



f[^^ = Ai / dx 

^0 [^ 



(39) 



The integral is only convergent if < 2K < 1. On the lattice convergence will be 

insured by a lattice cutoff ^ a as lower bound of integration. We can then make the 
substitution x = v arcsinh(w) / (27rT) and use partial integration to separate the cutoff 
dependent and independent parts. The cutoff independent part is expected to be the 
same for all K and given by 

p(i,con.) _ Ai^ r(l ~ 2K) [csc(7rX - i#) + csc(^A' + i#)] / 27rT y^-^ 

There is no bulk contribution in this order of perturbation theory. Taking derivatives 
with respect to h and T, respectively, the known results for the boundary susceptibility 
and specific heat are obtained [16, 17]. 

To obtain a correction to the bulk free energy we have to go to second order in Ai. 
As the bulk part is not influenced by the boundary conditions and we have already 
obtained the leading correction to the boundary part in first order we will use PBC 
to calculate this correction. In this case we have translational invariancc leading to 

/f ) =-^Jd^^ (cos + ^) COS (»^' 

We first consider the case T = where this correlation function can be obtained from 
(25) and (27). The imaginary-time integral is given by 

dr ^ K^T{2K-II2) 1 

(x2+w2r2)2^f vT{2K + l) a;4^-i ■ ^ ' 

This integral is always convergent because we assume x>c where c is a lattice cutoff 
of order a. The cutoff independent part of (41) is then given by 

For finite temperatures we can use again the conformal mapping of the plane onto a 
cylinder for the correlation function. The free energy then becomes 

/(2) = _M ('1L\'' r CO. '^dxdr 

^ \ V ) Jo Jo [sinh ^{x + wt) sinh ^{x - ivr)] ' 
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To evaluate this integral it is convenient to use the imaginary part of the retarded 
instead of the imaginary time correlation function and to introduce new variables 

Ui = X — vt, U2 = —X — vt leading to 

,,2 / rr.^ iK 



f (2) _ _ ^1 



sin(27rii') 



f 

Jo 



f) 

exp{iKhu/v) du 



(45) 



f 

Jo 



exp(— iA'/iw/w) du 



sinh^^ 

V 

This type of integral can be found to be given by [20] 

exp(iuz) 22^-1"' ' 



sinh^^ 



du 



sinh^^ ^1 



-B K-i- 



^-2K 



where B{x,y) 



ttT \ 2ttT' 
T{x)r{y)/T{x + y). This allows us to obtain the final result 



(46) 



(2)_ \l(2^T\ 

r(K- 



AK-2 



sin(27ri<')r2(l - 2K) 



(47) 



Kh 
27rT 



r {K 



V{l-K-i^)T{l-K + i^) ■ 
The other irrelevant terms in (36) with amplitudes a, /? can be expressed as 



2(9, 



=16 



(48) 



When we now shift the boson field according to (17), terms with an odd number 
of (/)-ficlds appear. Such terms do not contribute when calculating their expectation 
value in first order perturbation theory, because the free boson Hamiltonian consists 
of an even number of (?i-fields. The operators {dx(p)^,{dx(p)'^{dr(t>)'^,{dr(t>)'^ and 
(9^(/)7j)2, (9^(/)l)2, (9^c/)i?)((?^(/iL) independent of h will only yield a T^'-contribution 
to the bulk part of the free energy [14] and a T^-contribution to the boundary part. 
The T^-term in the bulk free energy is sub-leading compared to the T^-term in (35) 
and we will ignore it. The same is true for the T'^-tcrm in the boundary free energy 
which will be sub-leading compared to a T^-term which cannot be determined within 
the field theory approach as discussed at the end of subsection 2.3. 
We therefore have to consider only the following three terms 



AK'^h^ SAKh"^ 

where we have defined 

, a + 2(3 



167r2 



C 



2a- 12/3 



CKh^ 
2ttv^ 



The first term in (49) yields a temperature independent contribution given by 



JA,a 4^2„4 

The next term gives 



h" 



JA,b 



AT 3Kh^ 



J d^x{id,<i>f) . 



(49) 



(50) 



(51) 



(52) 
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The correlation function {{dxcff) at finite temperature can be obtained by taking 
derivatives of the correlation function in Eq. (23) and applying the conformal mapping. 
To evaluate (52) we will do a point splitting by a small parameter e in the correlation 
function. This leads to 



2v^L Jo isinh^i^e sinh2^(2a; + e) 



Doing the integral and then expanding in e yields 

(1) 3AKh^ 3AKh^ AKh^ 1 3AKh^ 

There are two different ways to interpret this result: One might think of e as a sort of 
lattice cutoff. Then one would say that the first term yields a correction to the constant 
in the bulk susceptibility which depends on this cutoff. However, we have already 
included all such corrections by using the exact value for the Luttinger parameter K, 
so this term has to be ignored. The second term shows that there is also a constant in 
the boundary susceptibility, however, it is also cutoff dependent. This means that we 
cannot determine this constant term within the field theory as expected from the more 
general discussion in subsection 2.3. The last two terms which are cutoff independent 
are therefore the only contributions we can really determine within the field theory 
approach. Another point of view is that we should use normal ordered operators. In 
this case we have to subtract the T — correlations from (53) thus exactly cancelling 
the first two terms in (54). In the limit e ^ only the two terms independent of e 
would then remain. 

The calculation for the third term in (49) is analogous and the result is 

(1) _ CKh^ CKh^ CKh^T^ 1 CKh^T 
Jc - 12^ '^LSttv ■ ^^^^ 

With the help of (50) we can summarise our results for the bulk part of the free 
energy as follows 

A7(2 ^ K^h^ ttT^ ^ r..Kh^T^ , , 

/bulk = eo - -— + (a + 2/3) — ^ - + (c, + 6/3) -—^ (56) 

XI sin(27ri^)r2(l - 2K)T {K ~ i^) T {K + i^) (2ttTV''-^ 



Av r(i-i^-i^)r(i-x + i^) V V J 

where eo is the ground state energy at zero magnetic field which is known from BA. 

The boTindary terms calculated by perturbation theory arc obtained for the model on 
the semi-infinite line, i.e., with one boundary. We therefore have to multiply these 
results by 2. The boundary free energy is then given by 

Fb = - Bh^ + + 0{T\T\ h^T\ ■ ■ •) (57) 



+ 



AiTT r(l - 2K) [csc(7r/s: - i-^) + csc(7rii: + i-^)] /27rr^ ^^"^ 



2ttTJ \ ' 2irT 



where Eq is the boundary ground state energy in zero magnetic field which wo will 
calculate by BA in the next section. The c;onstant B cannot be obtained within 
the field theory approach but will also be calculated in the next section by BA. 
0{T^, T**, h'^T^, ■ ■ •) denotes the regular terms which have been argued to be present in 
the boundary free energy based on the results for the free fermion model (30). These 
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terms remain undetermined here. The amphtudes Ai , a and /3, on the other hand, 
can be determined by comparing the field theory with the Bethe ansatz result as has 
been first shown by Lukyanov [14] for the periodic chain. We will show in section 4 
that the amplitudes obtained by comparing (56) and (57) with BA calculations for 
the open chain are consistent with Lukyanov's result. 



3. The Bethe ansatz solution 



In this section, leading and next- leading terms in a small-field expansion of the ground- 
state energy are calculated within the Bethe ansatz solution for the open chain. Using 
the Wiener-Hopf procedure, we obtain exact expressions for both the bulk and the 
boundary contributions. Section 3.2.2 contains an erratum to appendix B in [11]. 

In order to introduce notations and to make this article self-contained, the Bethe 
ansatz is shortly reviewed at the beginning, largely following [8]. We proceed to the 
thermodynamic limit afterwards. Leading and next-leading contributions to the bulk 
and boundary susceptibility in terms of the magnetic field are calculated subsequently. 



3.1. Coordinate Bethe ansatz for the open XX Z -chain 



The Hamiltonian (1) conserves the total spin in ^-direction = that 
the eigenvalue of is a good quantum number. Let us denote an eigenstate with 
= N/2 — M by |M). This eigenstate is a superposition of states with M-many 
J,-spins, which we denote by 

\ni,...,nM) = S~^...S-J0) 
where |0) is the fully polarised state with M = 0. Thus 

\M) = ^a(ni, . . .,nM)\ni . ..um) =■ ^a{n}\{n}). 

{n} {n} 

In terms of the coefficients a{n}, the eigenvalue equation H\M) = E\M) reads 



i ^ (a{n'} - Aa{n}) + ^A{N - l)a{n} = Ea{n}. 
{"'} 

Consider the simplified case 

1 < ni < . . . < hm < N, rij+i — nj > 1. 
Then the {n'} are given by the {n}, with = ria ± 1, so that 

M 

i> fc,no 



a{n} = exp 



and 

M 



^ a(ni, ...,na + l,.. .,nM) + a(ni, . . . ,71^ - 1, . . . ,nM) - 2Aa(ni, ...,na, 



= 2JEa{ni, . . . , um)- 
Inserting (60) yields 



M . 
E = J^cosfcj -I- JA ( -M-h 

3=1 ^ 



TV- 1 



(58) 
(59) 
(60) 

, njvf) 
(61) 

(62) 



To solve the full eigenvalue problem (58), two complications arise: 
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i) Flipped spins may be nearest neighbours. After acting with on a configuration 

with adjacent spins, amphtudcs arc obtained where the flipped spins are equal. 
These amplitudes must vanish. We thus first extend (59), 

1 < ni .. . < < iV. (63) 

and then solve the eigenvalue problem by requiring that (61) and therefore (62) 
still hold. This means that the following terms in (61) must be zero: 

a(. . . , n„ + 1, na + 1, . . .) - Aa(. . . , n„, n„ + 1, . . .); 
a(. . . , n„, ric, . . .) - Aa(. . . , n^, ric, + 1, . . .) 

which leads to 

a{. ..,na + l,na + 1, . . .) + a(. . .,na,na, • • •) - 2Aa{. . .,na,na + 1, . . .) = 0.(64) 

ii) We have to deal with open boundaries. Therefore, we extend the lattice to include 
the sites and + 1, where transitions to these sites are excluded: 

a(0,n2,...)- Aa(l,n2,...) =0 (65) 

a(ni,...,nM-i,A'+l)- Aa(ni,...,nM-i,A'') = (66) 

We expect that i) leads to a scattering phase in the amplitudes whereas ii) yields the 

"quantisation condition" of the ks. 

In order to solve i), Bethe [2] made the ansatz 

M 



PaTla 



(67) 



a{ni, . . . ,nM) = E ^(-P)exp 

PeTTM 

where the sum carries over all permutations of M integers. Then (64) reads 

^^A{P) ^e^'-'^p^+'^pi^+i)^ — 2Ae''^^<°+^' + 1^ e'(fepi"i+---+(fei'a+fej>(Q+i))no+---) _ 
p 

where Ua = Ua + I- Since the terms A{P) and A{PPa,a+i) have the same {n} 
dependence, (68) is fulfilled if 

A{P) ^e'('=^°+'=^(°+i)) - 2Ae"=P(«+i> + 1^ 

+A(PP„+i) (e'('=^(°+i)+'=^") - 2Ae^'=^« + l) = 0, 
from which we conclude that 

A(PP„+i) = -yl(P)e-^^^«-^("+i) (69) 
with the scattering phase 

ei(fej'a+fej'(c+i)) _ 2Ae'''^'(°+i) + 1 

g-lfpa,P(a + i) _ [ /YQ\ 

Qi(kp(a,+i)+kpo.) _ 2Ae''=^'<» +1 

From (69), one concludes that 

A(P) = c^S-o^^"'^". (71) 

If the chain was infinitely long, boundary conditions would not matter and one would 
choose the ks arbitrarily real (mod 27r). However, to carry out the thermodynamic 
limit properly and to obtain the boundary contribution to the ground-state properties, 
one applies boundary conditions for the finite system. Here, these are open boundary 
conditions (for periodic boundary conditions, cf. [8]). 
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Before continuing, let us parametrise the ka and 6ce,f3 by roots A^, such that (70) 
is fulfilled: 

pifca = sinh(Aa + i7/2) 

sinh(A„ - i7/2) ^ ^ 

eie(A..A,) ^ sinh (A„ - A,, + i^) 
sinh (Aa — A/3 — 17) 
In order to meet the requirement ii), the ansatz (67) is modified such that 
a(ni, n2, ...,nM)=J2 ^'(ei' • • • > ^m) J] ^"<^ «(A„-A,)+iE„ fep.n. (74) 

{e} P 

with ka = ea\ka\, Aq = ea|Aa| and we have defined the 2^ sets {e}, ej = ±1. Note 
that we have made use of (71) in (67). It is important to note that due to the sum 
over signs, (74) makes sense only for ka ^ 0, ±7r, that is A^ 7^ 0, ±m/2 for all a. 

The modified ansatz (74) does not influence (72), (73). However, it allows to 
meet the open boundary conditions (65), (66). Eq. (65) yields 



C(ei,...,ep„...,eM) ^ ei^^^-"^-") - Ae"=- 



From (70), 



+C(ei, . . . , -ep„ . . . , em) I 5^ ^-.oix,,+x,„) _ ^^-ife^x j = q. (75) 



such that (75) becomes 



C(ei, . . . , -epi, . . . , Em) 



= e5 E;3^„(e(Aa-As)+e(A<,+A^j)+2e(2A„)) 



C(ei, . . . , epi, . . . , Cm) 
with the unique solution 

C(ei,...,eM)=e-iS.<.3e(^«+A.3). 

We thus have the amplitude 

a{n} = ^^e5Sco(«(^c-A^)-e(Aa+A^))-iEa«(2Aa)+iEc,"afcPa_ 

p M 

This expression is substituted into Eq. (66) to obtain 

^iNkpM ^gifcpM g5 Z]„<M ^(Apa — Apm) — ^(Apq+Apm) 

^^-iTVfcpM ^g-ifepM _ Z]c<<M ^(Apa+ApM)-0(Apc-ApM) _Q jI'^y) 

Note that the fcs are parametrised by the As via (72). Using (76), Eq. (77) is written 
as 

a2(^+i)(Afc,7/2) a(2Afe,-7)_ (?M(Afc + i7)QM(-Afe - i7) 



a2(w+i)(Afc,-7/2) a(2Afc,7) gM(Afc - i7)(ZM(-Afc + i7) ^ 
with the definitions 

M 

a(A, ji) := sinh(A + i/x) , qni^) ■= Y\. ^iiih(A — Xj) . 

3 = i 



(78) 
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2N 



a(Afc,-7/2)_ 



a{2Xk, 7)a(Afc, 7r/2 - -f/2)a{Xk, 7r/2 - 7/2) 
a{2Xk,-j)a{Xh, -n/2 + j/2)a{Xk, -7r/2 + 7/2) 

^ i7)'?Af(-Afe - 17) 



(79) 



^AflAfe - i7)gM(-Afe + 17) ' 

which also foUows from the algebraic Bethc ansatz [10]. The Eqs. (78) are the above 
mentioned "quantisation conditions" on the As and therefore on the ks. 

The energy (62), including a magnetic field h along the S"^-direction, reads in 
terms of the As: 

M 



E =J 



sin^ 7 



N-1 



cosh(2Aj) — cos 7 
= N/2 - M. 



+ — r^cos7 -/iS^ (80) 

(81) 

The distribution of the roots has been investigated in [11]. There it has been 
shown that the equations (79) can be symmetrised by introducing the set of roots 
{vi,...,vn} := {— Ajv/2, • • • , — Ai, Ai, . . . , Ajv/2}, whose elements are distributed 
symmetrically on the real axis w.r.t. the origin. The vj are the N real solutions 
to the equations 

(82) 

Note that ±i7r/2, are solutions of (82). However, these solutions are not permitted 
within the Bethe ansatz as explained after Eq. (74). | This is a direct consequence of 
the broken translational in the open system. 



" a(Afe,7/2) ■ 


2N 


a{Vm,'!T/2 


-7/2)a(i;„,7/2) 


_ QNiVm 


+ 17) 


[a(Afc,-7/2)J 




_a{vm, -Tr/2 


+ 7/2) a(w„„ -7/2) 


qNiVm 


- 17) 



3.2. T = properties in the thermodynamic limit 

To pass to the thermodynamic limit, it is convenient to define the density of roots on 
the real axis, p{x). In [11], the following linear integral equation has been derived for 
p{xy. 



^{x, 7) + ^ Mx, 7) + ^x, TT - 7) + i?(a;, 27)] 

= p{x) + i?(a: - y, 2-y)p{y) dy, 

J-B 



where 



27ri'!?(a;,7) := 



2i sin 7 



d J sinh(a; + 17/2) 
dx sinh(a; — 17/2) 



(83) 



i-M.).(84) 



cosh 2x — cos 7 

This equation is valid including the order 0{1/N), that is, including the bulk and the 
boundary contributions. 

Equation (83) is a linear integral equation with two unknowns, B and p. In a 
first step, (83) is solved for i? = 00; in a second step, p{x) is obtained depending on 
the parameter B and the dependence of B on the magnetic field h is calculated. We 
will see that B = 00 corresponds to /i = 0, and a finite magnetic field h> induces a 

X Within the algebraic Bethe ansatz, the creation operators which are used to create eigenvectors 
by acting on a reference state are identically zero at spectral parameters 0, ±i7r/2. This can be seen 
directly from Sklyanin's work [10]. 
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finite B < oo. Finally, the susceptibility is deduced. This procedure is reviewed 
in [5]. 

It will be shown later that instead of dealing with p, all quantities we are interested 
in can be expressed more conveniently by g+{x) := 0{x)p{x + B). The calculation of 
these functions is done by Fourier transformation, 



1 

27r J_oo 



Let us first consider the case B = oo. It is straightforward to solve (83) in Fourier 
space, where 

sinn7rfc/2 

We denote the solution of (83) for i? = oo by po and find 

1 cosh7A;/4 cosh(7r/4 — 7/2)fc 



with 



2 cosh 7^/2' 27COsh7rx/7" 

Note that po(x)dx = 1/2 + l/(2iV), which corresponds to M = N/2 in Eq. (82). 
Putting this into Eq. (81), one has = 0. Thus a vanishing magnetic field 
corresponds to B = oo. This means that in order to perform a small field expansion 
\h\ <^ a (where a is some scale which is determined later), one has to expand 
asymptotically the energy at B ^ oo. 

To do so, we now consider the case B < oo, i.e., a finite magnetic field. Using 
(85) we can rewrite (83) as 

pix) = po{x)+ / k{x - y)p{y)dy (86) 
J\y\>B 

1 sinh(V2-7)fc 



"^""^ '~ 277 2cosh7fc/2 sinh(7r - ^)k/2 ^ ' " '^^^ (^^) 
We now introduce the functions 

p{x + B) =: g{x) = g+{x) + g-{x) 

g+{x) = e{x)g{x) , g-{x) = 9{-x)g{x) . (88) 
Then g{x) satisfies the equation 



(89) 



pOO pOO 

g{x) = po{x + B) + K{x-y)g+{y)dy+ k{x + y + 2B)g+{y) dy. 

Jo Jo 

The driving term po(x + B) can be expanded in powers of exp [B]. Since (89) is linear 
in g and po, we make the ansatz g = g^^^ + g^"^^ + . . ., where superscripts denote 
increasing powers of exp [B] . Then 

/■oo 

g^^\x) = [po{x + B)f^ + k{x- y)g^+\y) dy (90) 

Jo 

5(")(x) = [po(x + B)]("'+ / K{x + y + 2B)gX'-'\y)dy 

Jo 

+ K{x-y)g^^\y)dy. (91) 
Jo 
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Thus in each order, a hnear integral equation of Wiener-Hopf-type has to be solved. 
This technique is explained for example in [21, 22]. 
The first two orders of 5_|_(fc) read: 

tl\k) = G+{k) [po(A=)G_(fc)e-''=^]^'^ (92) 

fff (fc) = G+(fc) { [Mk) G.{k) e-"=^]f 

+ [Hik) g^'^ i-k) G- (k) e-2"^^] 1^'^ I , (93) 
where the indices ± arc defined by 

The functions G±{k) are obtained from the factorisation 1 — k = 1/(G+G_), they 
read 



- ln(7r/7 - 1) - ln(l - j/n) 

TT 



G (k) = V2(^^r(i-ifc/2) , 

+ r(l/2-i7fc/(27r))r(l-iA:(7r-7)/(27r)) ^ ' 

1 

" =2 

G_{k) = G+{-k). (96) 

We restrict ourselves to the calculation of which is sufficient for our purposes. 

The bracket [. . .J^^-* in (92) and the first term in (93) is evaluated using (94). We 
thus have to find the residues of 



1 cosh^ cosh(f - f) A; 



2 cosh ^' cosh ^ cosh 



(97) 



at the poles closest to the real axis in the lower half plane (the first term in the above 
line accounts for the bulk contribution, the second for the boundary contribution from 
(85)). This is done straightforwardly for the bulk part: the poles are located at 

/c^i) = -i(2n+l)7r/7, n = 0,l,.... (98) 

The relevant poles of the boundary part, however, depend on whether 7 < 7r/3 or 
7 > 7r/3: Poles at kn'' are found as well as poles at 

= -i2(2TO + l)7r/(7r - 7), m = 0, 1, . . . . (99) 

For 7 > 7r/3, the leading pole is k^^^, whereas for 7 < tt/S, the pole fcg^^ is leading. 
Finally double poles occur at 

- = A -i ' = + 2n + 3), m,n = 0, 1, . . . , (100) 

where the leading and next-leading ones occur for 7 = 7r/3.37r/5. In the ongoing, wo 
will concentrate on the single poles, i.e. 7 7^ 7r/3 in the leading order and 7 ^ 7r/3, 37r/5 
in the next-leading order. The extension of the results to the cases 7 = 7r/3, 37r/5 is 
discussed in section 5. We first consider the leading-order contributions (92), before 
proceeding to the next-leading corrections (93). The isotropic case 7 = is treated 
separately. 
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3.2.1. The leading orders By taking only the poles fcg^', kj^^ nearest to the real axis 



ao 



r/7 



-Tz/'yB 



2N 



«1 .g-7rB/7 



r/7 



+ 1 



-27rB/(7r-7) 



, // ^- (101) 

We could proceed analogously for 7 = tt/3 by evaluating the residue at the double 
pole [11]. However, it turns out to be more convenient to include the point 7 = 7r/3 
only at the very end, once the susceptibility has been calculated. The constants in 
(101) are given by 



ao = -G_(-i7r/7) 

7 

V2i . sm7rV(47) 

2i 

bi = tan7r7/(7r — 7) i27r/(7r — 7)). 

TT — 7 



cos(7r^/(47) - 7r/4) 



We can now compute := S^N and e := E/N from (80), (81): 



= 1/2- / p{x)dx + l/{2N) 

J-B 



e = -hs'' - 



J sin 7 



J 



'd{x,^)p{x)dx + — ( C0S7 + 



2 — cos 7 
TV 



We insert (86) into (103) to obtain 
s^ = ^^?+(0), 

TT — 7 



(102a) 
(1026) 
(102c) 

(103) 
(104) 

(105) 



which is an exact statement, including all orders 5'"'. It is convenient to calculate 
e — eo, where eo := e{h = 0) is the ground state energy at zero magnetic field. From 
(85), (104) one finds 



eo = — I C0S7 + 



= -cos7- 
4 Att 



J 
4 

J sin 7 
47r 



2 — cos 7 ^ J sin 7 
N 
J sin 7 



'&{x,^)po{x)dx 



f 



sinh(7r/2-7/2)fc 1 
2cosh7fc/2 sinh7rfc/2 A'' 



J(2 — cos 7) 



We use again (86) which yields 
4J7rsin7 



e - eo 



-hs 



cosh 7/0/4 cosh(7r/4 — 7/2)fc sinh (7r/2 — 7/2) k 
2cosh7fc/2 cosh(7r — 7)A;/4 sinh7rfc/2 

3+ (a;) 



dk 



(106) 



/ 

JO 



da; 



7 ^0 cosh(a; + -B)7r/7 
'^-(^i^)(0)+^f(0)) + «^^^^^^ 

ffi'^(3i7r/7)e-3-s/'r + O (e-'^-^/^gf^)^ 



^ (^i^)(0) +5^f (0)) + 5^ [(?i^^(W7) +?^^)(W7)) e--/- 



(107) 



where in the last equation we restrict ourselves to the given orders. Now B is treated 
as a variational parameter and is determined in such a way that 

^(e-eo) = 0. (108) 
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In this section we consider only the leading order in (107). Inserting (101), (105), 
(107) in (108), B is obtained as a function of h, 

B= -^In- (109) 

TT a 

27rJsin7(7r-7)G'+(i7r/'-,) 

"•=^ ^'''^ 

Thus a sets the scale for h (this scale a should not be confused with the amplitude a 
introduced in section 2). The restriction to the leading orders in exp[— _B] is equivalent 
to the leading orders in /i in the limit \h\ <C a. 

One now makes use of (110) to determine s^{h) from (105), and therefrom 
X{h) = ds^/dh. Inserting the expressions for G± from Eqs. (95), (96) we find 

Xbulk = 7 ^-r- — • 

(tt — 7)7r,7 sm7 

This result is well known, see, for example, [5]. The boundary contribution is given 
by 

- ^ sin7rV(47) 

' J(7r - 7)7r\/2 sin 7 cos(7r V (47) - 7r/4) 

(7r-7)2 7r-7Q r(l/2 + 7/(7r-7)) ^ ' ' ^ ' 
for 7 ^ n/3 with 

_ , ^sin7 r(l + 7r/(27)) 

a = 2J(. - j)V^—- (lllc) 

Note that the first term in (111&), which is independent of the magnetic field h, is the 
leading contribution for 7 > 7r/3 (pole closest to the real axis in (92)). For 7 < 7r/3 
the second term dominates. In how far the result (1116) also yields the next-leading 
corrections for the whole range of < 7 < tt cannot be answered at this point: so far, 
we have only focused on the leading poles. Next-leading corrections will be discussed 
in the next section. 

The constant contribution in (111&), as well as the pre- factor of the /i-dcpendent 
part, show divergences for certain values of 7. We will comment on these in section 5. 



3.2.2. The next-leading orders In this section, the Wiener-Hopf result (93) in (107), 
(105) is used to calculate next-leading corrections to (1116). This section also 
constitutes an erratum to appendix B in [11]. 

Let us start with the bulk contribution. We already commented on the calculation 
of the first term in (93) after that equation. The second term involves the residues of 
K{k) at the poles closest to the real axis. As can be seen from (87), poles occur at kn^ 
(cf. Eq. (98)) and at 

kP =-i^^l, 1 = 1,2,.... (112) 

TT — 7 

This situation is very similar to the leading order of the boundary contribution. 

The boundary contribution is treated analogously. The only difference is that the 
function g^_^\k) given in (101) contains two terms in 1/iV. The first term obviously 
yields exponents in analogy to the bulk. The second one has poles given in Eq. (99) 
which have to be combined with those in Eq. (112). 
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In the ongoing, we will focus on the amplitudes of the most important next- 
leading terms, that is, we consider the residues due to fcg"'^'' and feg^"* in the second term 
in (93). Calculating the residua mentioned above and putting everything together, 
one ends up with 

C-OA _|_ "'0,2 \ ^_3tB/7 _| ^0^3 ^_(7r/7+47r/(7i--7))B 

fi37r/7 k + 'm/^J fc + i27r/(7r — 7) 

"■1.1 _^ "'l.^ '\ g-37rB/7 _! "1.3 g-(7r/7+47r/(7r-7))B 

fc + i37r/7 k + in/jj fc + i27r/(7r — 7) 



1 

^2N 



+ 



k + i67r/(7r — 7) 



I ^1-3 \ g-67r/('r-7)S ^ ^1,2 ^-2(,r/7+7r/(7r-7))B 

fc + i27r/(7r — 7)/ A; + i7r/7 

i „ / .37r 
oo,i = — G-\ -1 — 
7 V 7 



ao,2 



1 TT^ /■ .TT 



tan ^ 



277r 27 V 7 

i 7r7 • \ ■ 27r 



^0,3 = —7 r tan G_ — i— ) G 



7r(7r + 7) TT — 7 



-1- 



_ . 2 sin7r/4 sin37rV(47) f .37r\ 

~ '7 cos(37rV(47) + 7r/4) V"' Y 7 

ai,2 = 7^ tan —G_ -1- 
27r 27 V 7 

°'i7 ^„„ 7r7 ^2 : 27r 



ai,3 = — ; — 7 tan G_ -i 

7r(7r + 7J TT — 7 \ TT — J 

, .2 377r / . 67r 
Oi 1 = 1 tan G_ — 1- 



TT — 7 TT — 7 \ TT — J 

,, = ^l(!^tanf G^_(-i^ 
7r(7r + 7) 27 V 7 

^1 7^7 ^2 f , 27r 



113) 
114) 
115) 
116) 
117) 
118) 
119) 
120) 
121) 
122) 



5i 3 = — tan G_ . 

At: TT — 7 \ TT — 7^ 

This expression for §2 is inserted into (107), where we now have to keep all the 
indicated terms. Then B as a function of h is derived. In section 3.2.1, we found 

that this relationship is the same both for the boundary and for the bulk in the 
leading order. This is no longer true when next-leading terms are considered. For the 
bulk we obtain 

„e-/^^.(l + A,(i)%A(^)'''""^') (123, 

ao,2 , G_(-i37r/7) tt - 7 00,3 

^1 = ^ I ■ I \ ' ^2 = —z ■ 

ao G_(-i7r/7) 27 ao 

Care has to be taken by inverting the relation h = h{B) in order to get B = B{h) 
for the boundary: this is done by performing a large-B-expansion in /i = h{B). The 
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next-leading terms in this expansion depend on the range of 7, whether 7 < 7r/3 or 
7 > 7r/3. In the final result B = B{h), the parameter 7 must enter uniquely. Then 
we find 

„e-/^ = .(l + ^.(^)%^.(^)""""^') (124) 

ai G-{—i'k/')) 27 ai 

Here terms of higher order, previously given in [11], have been discarded for the reason 
mentioned above. 

Combining these equations with (105), one finds 



'L,.('.) = ^^{G-(-»/7)i 



1+47/(77-7) 

tan-^G_ I -i- ) I -i^:— ) I - 1 (126) 



TT — 7 ^ 7r7 ^ f ■''^\ r^i2 f ■ '^''^ \ f h 



7r(7r + 7) TT — 7 \ 7/ \ it — ^ J \a 

27/(77-7) 



7 G+(i37r/7)\ //i' 

+ (27a,. + -a., + 7ai^^^j U 



T — 7 



These results correct those of [11]. Prom these expressions, Xbulk and Xb can be 
obtained: 



Xbulk(/i) = \ I ^ . {G_(-i7r/7)^ 
y 7r(7r - 7) t a 

3/1 7r^^o/.7r\ 2^ / .37r\\//i^^ 



tan^G3_ -i- +^G 



-1- 



Q:\7r 27 \ 7/ 3 \ 1 ) ) \oi. , 

(-i^) Gl f-i^) fi)""""'' (128, 

a7r(7r + 7j tt — 7 \ 7/ \ tt — 7/\a/ 



Xs(/i) = -ii 



1 fT^^Tfci 

27r(7r — 7) [ q; a \a J 



3 / 7 G+(i37r/7)\ / 

a \ 3 G+(i7r/7) J \a 



"-T f h\''-i ( h 



T-7 



5.^.5. Isotropic XXX-point The isotropic case 7 = (i.e. A = 1) is treated in the 
same manner as the anisotropic case 7 7^ 0. Therefore, we rescale (80) by Xj ^Xj. 
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This is equivalent to substituting fc — > A;/7 in Fourier space. Then 

1 



^^^^ "2coshfc/2 
PoW-W+2^2^(l + e-l^l/^). (130) 

Whereas the analyticity properties of the bulk contribution to (130) are qualitatively 
the same as in (85), the boundary contribution shows, besides poles, a cut along the 
imaginary axis. The functions g^' are determined as in the anisotropic case, where 
now the integrals in (92), (93) encircle the cuts of (130) and of 

e-|fe|/2 



. , , ,„e-''^^dfc (131) 
2coshfc/2 ^ ^ 

(that is the isotropic limit of (87)). Accordingly, the functions G± from Eq. (95), (96) 
now read 

°^W°^ r(i/2li./(2.)) --" 

_ _ 1 _ ln(27r) 
~ ~ 2Tr 27r 

G-{k) =G+{-k) 

Let us first consider the bulk contribution. We find 
5j'W+5?'W = G4«(^e--5^e- 

f §oe- 



kB2 ' fc ^ 

lM_:l+(3^lr^e--^+f3s^, k = 



(132) 



with 



* = ^Gi,0) ft = 

Inserting (132) into (107) and performing the maximisation condition (108), one 
obtains 

e--- = - f 1 + + "2^^^^^ + (133) 

a\ Inh/a \n^h/a \n^h/a) ^ ' 

with 

.Att^ 1 

ai = - 1 = -- 

oo 4 

./327r3 1 

oil = 1 = --^ 

ao 8 

TT^ 1 

as = — (i/337r + i/?i - 2/3o - i7r/J2 Iutt) = - - 
ao » 

Combining (132), (133) and (107), one obtains for the bulk susceptibility 

„x 2ai „ ln|ln/i/a| \ 
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with 0:4 = 5/16, Xbuik(^ = 0) = J/tt^, and a = \/2n^J^. We now set a = /ioi5, 
where b is determined such that the term ~ ln~^/i/Q: in (134) is absorbed into the 
1/ ln(/i,//i,o)-term. This prescription fixes the scale uniquely. One finds (5 = exp(— 5/8) 
and 

/lo = \/2^e-9/^ (135) 
Then Eq. (134) reads 

_ 1 / 1 ln|lnV/io| 

Xbuik-j^(^ - 21nV/Jo " 41n2V/io 
Let us now consider the boundary contribution. In (92), the [. . .]_|_-bracket yields 
contributions 0(cxp [—const. _B]) from the poles, and algebraic contributions duo to 
the cut. The exponential contributions are clearly sub-leading in comparison to the 
algebraic ones, so only the latter are calculated in the following. Using the expressions 
for (j±(fc), we find (omitting the bulk contribution) 
G+(0)(ai/B + a2(ln-B)/S2 + 03/52) 

+0((lnB)/B3,l/B3), fc = o 

iaiG+(fc)/(fcB2), fcT^O 

1 \/2 1 / 1 , , 

= 71^' = "i^' = 71^ (^ln2- -ln(27r) 

Leading and next-leading contributions are here already contained in g^^^ , so that we 

do not consider further corrections stemming from 'g^^ . Prom (105), (107), (108) we 
obtain 

B = - - In - (136) 

These equations are obtained by those from the anisotropic case, (109), (110), by 
scaling B jB and sending 7^0 afterwards. Carrying out the same steps which 
lead to (111&), one finds the boundary contribution 

sUh) = -U ^ + i^^i^") + o(ln- h) (138) 
if 1 ln|lnV^o| 1 \ / 1 \ , , 

ho = a/\/2 = JiT^/TT/e . (140) 

The scale ho has been chosen such that in (138), no terms 0(ln~^ h) occur. Note that 

ho ha. The results (139), (140) agree with the TBA-work by Frahm ct al. [23] for 
T = 0. Furthermore, agreement is found with [16, 17, 24], where scales which difi'er 
from ours (140) by a constant factor were used. 



4. Determining the amplitudes by comparing field theory and Bethe 
ansatz 

Apart from the amplitudes Ai,a,/3, also the bulk and boundary part of the ground 
state energy as well as the constant term in the boundary susceptibility have been 
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left unknown in our final field-theory result (56) and (57). These quantities have been 
calculated in the last section by Bethe ansatz. In terms of the Luttinger parameter K 
we can write the ground state energy as 

J J r sinf sinh^l^ ^^^^^ 

e^""- = - _ cos — - — / dk r#^TT^ (141) 

4 K sinhf cosh(^=^ ^ ' 



and 



i?o^ = ^f2 + cos^) (142) 



K 



— sm — smh— — / dk— 
Stt K 2KJ_^ si 



cosh cosh ^j^' 



sinh ^ cosh cosh ^^^^^ 



The constant B in Eq. (57) can also be expressed in term of the Luttinger parameter 

B=-^J_!!^. (143) 
2Trv 2V2 cos 

Next we determine the amplitudes. By taking the limit T — > of Eq. (40) we 
obtain 

/ , ^ 2K-1 

^(i,conv) ^ _x^K^K^{-2K) sin(i^7r) ( ^ J . (144) 

The amplitude Ai of the umklapp term can now be found by comparing this with the 
BA result in Eq. (1116). This leads to 

2K-2 

KT{K)suvtt/K r(l+2^' 



Ai 



7rr(2 - K) 



2v^r 1 



K 



2K-2 



(145) 



and agrees with Lukyanov's result in Eq. (2.24) of Ref. [14] where our amplitude Ai is 
related to his amplitude A by Ai = v\/2'!t. Here we have determined Ai by comparing 
with a boimdary quantity. Wc can. of course, also compare the zero temperature limit 
of the bulk contribution (47) with (126). This yields the same result as it should be. 

The amplitudes of the terms with scaling dimension 4 represent a more subtle 
problem. From the zero temperature Bethe ansatz solution (126) we can only obtain 
the amplitude of the /i^-term in the free energy. According to (56), however, this 
amplitude depends on a and (3. One possibility would be to calculate the amplitude 
of the /i^T^-term in /tuik from a finite temperature BA solution. Although such a 
solution exists, an analytic formula for this term is difficult to obtain. We therefore 
follow Lukyanov's idea to obtain both amplitudes from the ft,^-term in /b„ik alone, by 
using symmetry arguments. The coefficient of the /i^-contribution in (126) consists 
of two terms. The first of these terms vanishes at the free ferrnion point due to the 
tan- function. § We can identify this term as the one determining the amplitude a, 
because a is associated with a mixing of left- and right-movers and therefore has to 
vanish at the free fermion point. The other term then determines /3 leading to 

a = -TT w tan ^ (146) 

§ This does not come as a surprise: The term has been obtained from the second bracket in (93). 
However, at the free fermion point 7 = 7r/2, the integration kernel k vanishes identically, k\^—t^/2 = 0- 
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( 3A- \ p3 ( _J_\ 



and 



2A'-2 J \2K-2 

These results also agree with Lukyanov's formula (2.24) in [14] with a = —27t^vX+ 
and /3 = -2'k^vX-. 

5. Divergent amplitudes and logarithmic corrections 

From (56) and (57) the bulk and boundary susceptibilities for zero magnetic field can 
be derived by taking derivatives with respect to h leading to 

,_(T, = 0) = ^ - + 6/?)^ + ^r^(l/2 - K) 

xr2(l + K) sm{2nK) [*'(1 - K) - -^'{K)] ( ^ J (148) 

and 

KT KT(K)T(3-2K)(tt^ ~2m'(K)) f2TrT\'^^~l , 

X.(T, . = 0) = 2S - + A. ^J.(Ll/i.)r(2-i.) (— ) ^'''^ 

Interestingly, the constant term 2B is closely related to the constant in the bulk 

susceptibility (see (143)). The two terms become equivalent at the free fermion point 
as expected from the calculations in subsection 2.3. On the other hand, for T = 0, 
bulk and boundary susceptibilities are given in Eqs. (128), (129). 

For certain anisotropies, terms in x(T = 0, h) and x(r, h = 0) show divergences. 
In this section, we will show that these divergences cancel and give rise to logarithmic 
corrections. In the first part, we will focus on x(r = 0,h), whereas in the second 
part, x(r, ft. = 0) is treated. In order to keep contact with the notations introduced 
in sections 2 and 3, we will write our results in terms of 7 = arccos A for the T = 
case, whereas K = -k/{-k — ^) will be employed for /i = at finite T. 

5.1. Logarithmic corrections in x{T = 0,h) 

Let us come back to the enumeration of poles encountered within the Wiener-Hopf 
procedure in the leading and next-leading orders, Eqs. (98,99,112). Combining the 

relation B = B{h), Eq. (109), with the functional dependence of the energy on g+ , 
Eq. (107), we can read off the exponents of /i-dependent contributions. Namely, the 
poles (98) lead to terms 

n = 0,l,... (150) 

both in the boundary and the bulk susceptibility. Furthermore, the km^ from (99) 
lead to contributions 

~ h — — — , m = 0, 1, . . . (151) 

in the boundary part. As far as the next-leading order is concerned, one determines 
the exponents to the bulk susceptibility due to the second term in (93) from the 
combination of poles (98), (112). This results in terms 

^ f^2in+2m+l) ^ ^2(n+2^l) ^ n, m = 0, 1, . . . ; Z = 1, . . . . (152) 
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The next-leading order of the boundary susceptibility shows, apart from terms similar 
to (152), additional contributions of the form 

^ f,^n+l+^(2m+l) ^ ^^(2^+1+20-1^ n, m = 0, 1, . . . Z = 1, . . . . (153) 

which result form combining (99) with (98) and (112). Note that the pre-factors of 
the leading terms in Eqs. (150,151,152,153) are those which have been determined in 
section 3. By carrying out the Wiener-Hopf procedure to higher orders, it would be 
possible to obtain further pre-factors, at the expense of more and more cumbersome 
calculations. 

The enumeration of possible exponents in (150,151,152,153) is far from complete: 
These exponents have been found from extrapolating the first and second order results 
of the perturbation expansion of 5+, while using the first-order result (109) and the 
leading orders in (107). In higher orders of the 5f+-expansion, these exponents mix 
with each other and further orders appear. 

So far, we excluded the double poles in Eq. (100), which occur when two poles 
of the types (98,99) coincide, that is, when exponents stemming from (150,151) cross. 
Again combining Eqs. (109), (107) we obtain logarithmic contributions 

~/i^"ln/i, n = 0,l,... (154) 

to the boundary susceptibility. The same happens if exponents in the next-leading 
order contributions cross, therefore yielding logarithmic contributions of the form 
(154), with n > 1 there, for both the bulk and the boundary. In order to obtain the 
coefficient of the logarithmic terms, it suffices to know the amplitudes of the terms 
whose exponents cross: They show divergences at the crossover points, which cancel 
to yield the logarithms. 

Consider the constant term in XB{h), cf. Eq. (129). It displays poles at 

^" = 4;^T3 ^'''^ 
with an accumulation point at 7 = 0, whereas the coefficient of the /i-dependent 
contribution with exponent 27/(77 — 7) — 1 has poles at 

2m +1 

with an accumulation point at 7 = tt. Obviously, 7„ = 7^ for n = = m, which 
means 70 = n/3 (the case n = m = —1, 7_i = n, is excluded here). By setting 
7 = 7r/3 + e and expanding in e, the two divergences cancel and one obtains 

1 //ie^-''/2\ 
XB{h)U=./s=-^,ln[^^^j. (157) 

As far as the other poles (155), (156) are concerned, note that they coincide with 
(100) for m = or n = there, respectively. This means that we expect higher-order 
terms, whose coefficients show poles at (155), such that logarithms of the form (154) 
occur, with n > 1 there. 

A similar crossover at 7 = 7r/3 happens in the bulk part in the next-leading order, 
when the exponents 2,47/(77 — 7) cross in Eq. (128). Again divergences in the two 
corresponding amplitudes cancel, as can be seen by expanding 7 = -I -|- e. Then one 
finds 

1 2 /ie~^/^ 

Xbullcl.=./3 W = ^ - (158) 
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Finally, let us consider the additional crossover in the next-leading order of Xsih)- 
In Eq. (129), three exponents have been identified: 0, 2"//{t: — 7) — 1, 2. Contrary to 
the bulk case, there are now two crossover points: 7 = n/3 and 7 = Stt/S. At 7 = 7r/3, 
the (leading) logarithmic contribution has been identified in Eq. (157). At 7 = Btt/S, 
a (next-leading) logarithmic contribution appears, which is 



XBih)=Xo \ ^ + h 



2 , 135 (3 + \/5)r3(4/3) 

M 7r3/2(3 + V3){5 + 2\/5)r3(ll/6) 

(5-|-y5)7rr2(ll/6)A 



-181n/i-31 + 391n2-|-91n 



r2(i/3) J 



with 



405 (2 + ^/3)(3 + ^/5)^3(4/3) _ 27 JlW3)(3W5^\ \ .^^g. 
64 (3 + V3)7ri/2(5 + 2V5)r3(ll/6) 10 7r(3 + V3)(5 + 2^5) J J 

^0 = 1^^^(5 + 2v/5)(5 + ^)-^/2. (160) 

This completes our discussion of the leading and next-leading terms in the 
susceptibilities on the basis of the Bethe-ansatz solution at T = for the anisotropic 
XXZ-case. 



5.2. Logarithmic corrections in x(T, /i = 0) 

Let us start our discussion for finite temperatures with the bulk susceptibility. The 
T^-term in Xbuik contains the amplitude a which is defined in (146). a is divergent for 
K = {2n+l)/2n with n = 1, 2, • • -.H Obviously Xbuik does not diverge at these points so 
this divergence has to be cancelled by other divergent terms. In fact, the T'*^~'*-term 
has also divergences a.t K — (2m + l)/2 with m = 1, 2, • • • due to the r'^{l/2 — K)- 
factor. The only point where both terms simultaneously diverge is therefore if = 3/2 
and we will show that these divergences indeed cancel each other here and lead to 
a logarithmic correction. But what happens to the other divergences? Let's first 
consider the divergences in the r''^"4-term. ki K ^ (2m + l)/2 the temperature 
dependence becomes T''™-^ = r^, T^", • • •. Terms like /i^T^", • • • appear in 

the free energy in perturbation theory in the irrelevant terms with integer scaling 
dimension. For example, in n-th order perturbation theory in the {dx(j)Y operator 
we will get a h^T^'^~^-tei:v[\. These higher order corrections will interfere with the 
T^^^*-term and cancel the divergences -aI K — (2m + l)/2. The result will be a 
rp4m-2 Y^j^ behaviour at these points. As the T^-term dominates for K > 2 these are 
just logarithmic corrections in next leading contributions. 

Much more interesting are the divergences at K = {2n + l)/2n in the T^-term. 
They yield In T-terms and the points where this happens become dense for K 1 
(A ^ 1). The divergence at if = 3/2 in the T^-term is cancelled by the divergence 
in the T^'^~'*-term. In general, the divergence at K — {2n + l)/2n is cancelled by a 
j'4nif-4n_^gj.jjj which One obtains in 2n-th order perturbation theory in the cos(2^/i?)- 
operator. 

II The amplitude /3 does not show any divergences. 
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The only case where we can give an explicit result for the logarithmic term is 
K = 3/2. Here we find 

1 142 + 247+ 121n(3888)- 21C(3) 2 



Xbuik(r, /i = 0) = 



243\/37i 



-T^ (161) 



where 7 is Euler's constant and ({x) Riemann's zeta function. 

In the boundary susceptibility similar things happen. Here even the constant term 
(143) as determined by BA shows divergences at K = (4n— l)/(4n — 2) = 3/2, 7/6, • • • 
due to the cos-term in the denominator, see also Eq. (155). The T^^~'^-term, on 
the other hand, has a divergent amplitude for K = {2n + l)/2. Again all these 
divergences have to cancel each other leading to logarithmic corrections. For if = 3/2 
we can derive an explicit result 



XBiT, h^0)= -^ + ^^ + ^- + ^H^m-7a3) JnT ^ ^^^^^^^ 



(162) 



4V37r2 V37r2 

Here the TlnT-term stems from the term linear in T whose amplitude also diverges 
at K = 3/2. 

We summarise our results in Fig. 5.2. 
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Figure 1. Crossover of critical exponents, depending on the anisotropy A, in the 
leading orders. Black lines are those exponents which occur both in the bulk and 
boundary susceptibilities (namely and 2). The leading non-integer exponent in 
the boundary contribution is dashed blue {2K—3), the one in the bulk contribution 
dotted blue {4K — 4). Green is the exponent 1 , which only occurs at finite T in the 
boundary contribution. All crossover points are denoted by circles. The crossover 
between 1 and 2K — 3 at A = has no consequences, since the corresponding 
amplitudes vanish. At all other crossover points, logarithmic corrections appear. 
The vertical lines denote A = (1 - \/5)/4 (7 = 37r/5), A = and A = 0.5. 
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6. Comparison with numerical results 

To check our analytical formulas, in particular those for the boundary susceptibility, 
we use numerical data. For T = we have solved Eq. (83) numerically (for more 
details see [11]). This allows us to obtain XB{h) and to compare with the analytical 
formula (129) if \h\ <^ a. Such a comparison has already been performed in [11] in the 
repulsive regime < A < 1 verifying our analytical result for the constant term and for 
the term which involves a fractional power of h. These are the dominant contributions 
for 1 < iiT < 5/2. For K > 5/2 the leading term is quadratic in h. To check also this 
term we present in Fig. 2 numerical data in the attractive regime. Another point worth 
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Figure 2. a) Boundary susceptibility XB{h). Comparison between a numerical 
solution of the BA equations (black solid lines) and the analytical solution 
for h/J <^ a in the attractive regime (blue dotted lines), b) The difference 
XB{h) — X-b(O) on a double logarithmic scale. For A < cos(37r/5) —0.31, the 
leading h-term has constant exponent 2, whereas for A > cos(37r/5), the exponent 
depends on 7, according to Eq. (129). 

checking are the two anisotropics K — 3/2 and K = 5/2 where we have obtained the 
explicit formulas (157) and (159) for the logarithmic corrections. These are compared 
with numerical data in Fig. 3 and Fig. 4, respectively. Next, we turn to finite 
temperatures. The bulk quantities can be calculated numerically based on the Bethe 
ansatz solution in the quantum-transfer- matrix approach [7]. A comparison between 
these data and Eq. (148) in the repulsive and in the attractive regime is shown in 
Fig. 5 and Fig. 6, respectively. Figs. 5b) and 6b) show xiT, h = 0) — x{T — 0,h ^ 0) 
in double logarithmic plots. Excellent agreement is found for both the amplitudes 
and the exponents of the leading T-corrections, thus confirming (148) first derived by 
Lukyanov [14]. In particular, note the dependence of the exponent of the leading T- 
correction on the anisotropy in the repulsive regime. Fig. 5b), whereas this exponent 
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Figure 3. Numerical data for the boundary susceptibility (black solid line) at 
K = 3/2 (i.e. 7 = tt/3) compared to formula (157) (blue dotted line). The inset 
shows XbW — clnh, where c = —l/\/3n^. According to Eq. (157), this yields 
the scale of the logarithmic divergence. 
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Figure 4. Numerical data for the boundary susceptibility (black solid line) at 
K = 5/2 (i.e. 7 = 37r/5) compared to formula (157) (blue dotted line). For the 
plot in the inset, the constant Xsih = 0) has been subtracted. Note the double 
logarithmic scale in the inset. 
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Figure 5. Comparison between numerical data (black solid lines) for the bulk 
susceptibility in the repulsive regime and formula (148) (blue dotted lines) for 
7 = 0.5,0.8, 7r/3, 1.3, corresponding to A = 0.88,0.7,1/2,0.27 (from bottom to 
top in a) and from top to bottom in b), c)). In b), the constant x(T = 0, ft. = 0) 
has been subtracted; the inset c) is a zoom of b). 

For the boundary susceptibihty at finite temperature no Bethe ansatz solution 
is known today. We therefore have calculated this quantity by using the density- 
matrix renormalisation group applied to transfer matrices (TMRG). This method is 
particularly suited because the thermodynamic limit can be performed exactly and 
no finite size corrections disguise the boundary contributions we are looking for. This 
method has already been used in [11] to obtain data in the repulsive regime and 
the reader is referred to this article for more details about the TMRG. In Fig. 7 
we show data for the boundary susceptibility in the attractive regime. Interestingly, 
Xb first decreases when starting from infinite temperature down to a temperature 
Tg which depends on anisotropy before it starts increasing at lower temperatures. 
This behaviour can be understood in terms of a crossover from antiferromagnetic to 
ferromagnetic short-range order in the bulk with increasing temperature, as studied in 
[25]. According to [25], a change in the dispersion relation of elementary excitations 
from linear to quadratic behaviour occurs when increasing the temperature from 
T < To (7) to T > To (7). Note that the critical temperature To (7) depends on the 
anisotropy. 

We expect that the crossover temperature between these two regimes is described 
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Figure 6. Comparison between numerical data (black solid lines) for the bulk 
susceptibility in the attractive regime and formula (148) (blue dotted lines) for 
7 = TT - 0.5, TT - 0.8, 27r/3, TT- 1.3, corresponding to A = -0.88, -0.7, -1/2, -0.27 
(from top to bottom in a) and b)). In b), the constant x(T = 0, h = 0) has been 
subtracted. 
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Figure 7. TMRG data for xb (T) in the attractive regime. The black diamonds 
represent the zero temperature result known from BA. The black dots denote the 
minima To of xb{T). In the inset we show that these minima (black dots) are 
well described by formula (163) (red dots) with A = 1.245 obtained from a fit. 
The dashed lines are a guide to the eye. 
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by the same formula for TqIt) a-s for the bulk, which reads [25] 

„ , , ,sm7 TrfTT — 7) 

= A -tan \ 163 

7 27 

Here A is determined by a matrix element which is different for boundary and bulk 
susceptibility. This matrix element cannot be calculated by BA. In the inset of Fig. 7 
we show that this formula indeed describes very well the minima which occur in xb 
where A is obtained from a fit. 

Finally, we have a closer look at the low-temperature regime of Fig. 7 and compare 
the TMRG data with (149). Fig. 8 shows excellent agreement confirming that xb 
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Figure 8. Comparison between the TMRG data (circles) and (149) (solid lines) 
for XBiT) in the attractive regime. The black diamonds denote the T = BA 
result. 

depends linearly on temperature in the attractive regime for T ^ 1. 
7. Conclusions 

We have studied low-energy thermodynamic and ground-state properties of the open 
spin-1/2 XX Z-cham by combining the bosonisation technique with the exact Bethe 
ansatz solution in the critical regime — 1 < A < 1. Bosonisation has allowed us to 
obtain a low-temperature, low-field expansion of the free energy in terms of T, h for 
both the bulk and the boundary parts of the free energy with unknown coefficients. 
We have argued that the coefficients of the regular terms in the boundary part (these 
are those powers of /i, T which appear also in the bulk part) involve the short-distance 
lattice cutoff and therefore cannot be obtained from field theory. All leading terms in 
the bulk part, on the other hand, are cutoff independent and their coefficients could be 
determined from Bethe ansatz by expanding the ground-state energy in terms of the 
magnetic field h. For the boundary free energy only the coefficients of two terms turned 
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out to be cutoff independent. f These are a h'^T- and a /i^T^^~^-term. Fortunately, 

these two terms combined with a ft,^-tcrm obtained by BA give a complete low-energy 
description of the boundary susceptibility xb for all anisotropics. 

We used the Wiener-Hopf procedure within the Bethe ansatz in the 
thermodynamic limit to obtain the ground-state energy and a systematic low-field 
expansion of the susceptibility at zero temperature. The coefficients of the leading 
and next-leading terms in this expansion have been calculated explicitly for both the 
bulk and the boundary contributions. The possible exponents of higher order terms 
have been classified in a systematic way. 

We compared our results with numerical data obtained by various techniques: The 
T = Bethe ansatz equations involved solving linear integral equations numerically, 
the finite-T-behaviour of the bulk susceptibility was calculated in the quantum transfer 
matrix approach, and the finite-T-behaviour of the boundary susceptibility was 
obtained from the density-matrix rcnormalisation group applied to transfer matrices. 
In all cases, we found excellent agreement between analytical and numerical results. 

We identified several crossover phenomena. Depending on the anisotropy, the 
crossover of scaling dimensions of the irrelevant operators in the low-energy effective 
Hamiltonian leads to a crossover of critical exponents of T, h in the free energy. 
Associated with this crossover are divergences in the corresponding amplitudes, which 
cancel to yield contributions logarithmic in T, h at the crossover points. Whereas this 
happens in the next-leading orders in the bulk, this affects the boundary contributions 
both in leading and in next-leading orders. Furthermore, we observed a crossover in the 
boundary susceptibility in the attractive regime, due to competing antiferromagnetic 
or ferromagnetic ordering tendencies at low and high temperatures, respectively. 

We think that the rich results of the boundary behaviour analysed in this article 
stimulate further research in direction of an exact treatment of the boundary free 
energy at finite temperatures. 
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